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1. Introduction 

We present a first-principles calculation of off-forward matrix elements of operators that mea- 
sure the vector and scalar form factors of the pion, utilizing chirally improved fermions ||]. The 
vector form factor F n is defined by 

(K + ( Pf )\V^ |7T + (p,-))cont = (pf + PiinF^Q 2 ) , (1.1) 

where Q 2 = {pf — pij 1 is the space-like invariant momentum transfer squared, and = 2/3 uy^ u — 
1/3 dy^d is the vector current. Electromagnetic gauge invariance implies the constraint F n (§) = 1. 
The mean square charge radius is defined as 



The scalar form factor T n is given by the matrix element of the scalar operator 

(7l + (p f )\m u uu + m d dd\7l + (p i ))=r jl {Q 2 ) . (1.3) 

From chiral perturbation theory one knows that the scalar form factor at Q 2 = is the so-called 
sigma term which behaves like r^(0) ~ M 2 near zero momentum transfer. The scalar radius (r 2 )^ 
can be obtained from 

TM- l -\^ + em. 0.4) 

For a detailed discussion of chiral perturbation theory in this context see [BQ. 



s • (1-2) 

Q 2 =0 



2. Strategy 

In order to compute the form factors, we need to evaluate off-forward matrix elements (at 
several transferred momenta), i.e., the expectation values of 

Tr( £ e^S(0,y)O £ e-**S(y,x) y 5 S(x,0) y 5 ) , (2.1) 

where S(y,x) is the quark propagator from x to y and O denotes the operator inserted at y. We use 
the notation p = pj, r = p, and momentum transfer q = p — r. These matrix elements have been 
evaluated by using the sequential source method where the sequential propagator can be easily 
computed by an additional inversion of the Dirac operator for each choice of the final momentum. 
We extract physical matrix elements by computing ratios of 3-point and 2-point correlators: 

R(t z . Da] _ ( P (t-,P)O(r,q)P(0-,r)) / (Pft ;p)f (0;p)) (f (r;p)f (0;p)) (P(t - T;r)f (0;r)) 

U ' P ' q; <P(f;p)P(0;p)) V <^;r)P(0;r))(P(T;r)P(0;r))(P(f-T;p)P(0;p)) ' 

(2.2) 

where P denotes the pseudoscalar interpolator. We keep the sink fixed at time t and change the 
timeslice z where the operator O sits (scanning a range of timeslices). R(t,z;p,q) should exhibit 
two plateaus in z, for <C z -C t and t <^z <T. Note that the terms under the square root become 
trivial for q = (0,0,0). 
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16 3 


x32 


0.74 GeV 


1.04 GeV 


1.28 GeV 


1.48 GeV 


12 3 


x24 


0.98 GeV 


1.39 GeV 


1.71 GeV 


1.97 GeV 



Table 1: Values of the nonzero momentum transfers. 



We choose for initial and final momenta |py| = |p,| such that Ef = Ej and the transferred 4- 
momentum is given by Q 2 = |q| 2 . In this way one achieves for the electromagnetic form factor [|[] 
the cancellation of the kinematical factors in ( |1.1[ ) and 

(7r(p/)|V M |7r(p,-)}iatt = 7 == (ft(p / )|y y |7z;(p0)cont . (2.3) 

Indeed, with this choice of momenta, and using the pL = 4 component of the vector current, the 
overall factor becomes „ „ 

Ef + E: 

\ 1 = 1. (2.4) 

2yJ EfE t 

This cancellation is unfortunately no longer possible in the calculation of the scalar form factor, 
where a remaining multiplication by the quantity 2 ^jEfEj is still needed. 

We have made the choice |pj| = |p ( | = V2po, where po =2n/(aL) denotes the smallest 
(spatial) momentum unit. Fixing the final (sink) momentum to p/ = (l,l,0)po> we g et twelve 
values for the initial momentum. They give rise to four nonzero (and equidistant) values for the 
momentum transfer, Q 2 = 2np\ for n = 0, 1, 2, 3, 4 (see Table [[]). We could also consider larger 
values of the module of the initial momentum. In this case, however, the statistical errors become 



much larger, and moreover we could not use Eq. (2.4) anymore. A change of the momentum at the 



sink would require the expensive computation of a new set of sequential propagators. 

In the Monte Carlo simulations additional numerical noise due to the momentum projection 
of the sink and of the operator onto nonzero momentum transfer q arises. The 2-point correlators 
may even become negative (within the statistical errors) on timeslices near the symmetry point. A 
reasonable choice was to put the sink at a timeslice smaller than T/2, such that the 2-point functions 
have smaller errors. In particular, we have put the sink at timeslice t = 8, while the source sits at 
timeslice t = 1. 



3. Simulation parameters 

We use the chirally improved Dirac operator |Jlj, g], which constitutes an approximate 
Ginsparg-Wilson operator. The gauge action is the Luscher-Weisz tadpole improved action at 
/3 =7.9, corresponding to a lattice spacing of 0.148 fm (determined from the Sommer parame- 
ter). Although we expect improved scaling properties, here we do not study scaling behavior. The 
hadron interpolators are constructed from Jacobi-smeared quark sources and sinks. 

We have done simulations using two volumes, 12 3 x 24 (200 configurations for each mass) and 
16 3 x 32 (100 configurations for each mass). This corresponds to spatial lattice sizes of ah ~ 1 .8 fm 
and ah ~ 2.4 fm respectively. 

We work at several values of the bare quark mass; the corresponding meson masses are given 
in Table ^[ The jackknife method is used for estimating the statistical errors of our results. 
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16 3 x 32 


16 3 x 32 


16 3 x 32 


12 3 x 24 


12 3 x 24 


m 


0.02 a' 1 


0.04 a' 1 


0.06 a' 1 


0.04 a' 1 


0.06 a- 1 


m n 


342 MeV 


471 MeV 


571 MeV 


414 MeV 


575 MeV 


trip 


845 MeV 


895 MeV 


941 MeV 


952 MeV 


976 MeV 



4. Results 



Table 2: Pion and rho masses at the various quark masses simulated | 



In Fig. [j]a we show an example for the observed plateau behavior of the ratio (12). Fitting the 
central points of the plateau then provides values for the electromagnetic form factor of the pion. 

When comparing the results to physical, renormalized quantities we have to multiply with 
renormalization factors relating the raw results to the MS-scheme. For the chirally improved Dirac 
operator these have been determined in [^] and all results we show are already converted to the 
MS-scheme. Although the vector current is pointlike and not conserved, the value Zy = 0.9586(2) 
turns out to be close to unity. The resulting F n (0) is not constrained to unit value, but comes very 
close as may be seen, e.g., in Fig. [T]a. This fact is also obvious from Fig. |I]b, where we plot the 
inverse of the vector form factor determined for all the momentum transfer values and quark masses 
studied on the 16 3 x 32 lattice. 

The electromagnetic pion form factor in the time-like region is dominated by the p -meson. 
In the space-like region it may therefore be well approximated by a monopole form, F % (Q ) = 
m v I (m v + <2 2 ), such that 1/F n has the linear behavior indeed observed in Fig. [l]b. In that approxi- 
mation the mean square radius is 6/m v . 

As the vector meson dominance (VMD) model is just an approximation one expects correc- 
tions due to other resonances and more-particle channels in the time-like region. These may be 
taken into account by additional pole terms. The data (cf. Fig. [j]b) does not really require such a 
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Figure 1: (a) L.h.s.: The ratio R(t , T;p,q) leading to the electromagnetic vector form factor at zero momen- 
tum transfer, for the 16 3 x 32 lattice at m n = 470 MeV. (b) R.h.s.: The inverse vector form factor for the 
three available masses on the 16 3 x 32 lattice. For presentation purposes the data points for different masses 
are slightly shifted horizontally. 
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Figure 2: The mean square radius for the three available masses on the 16 3 x 32 lattice for (a) the vector 
(l.h.s.) and (b) the scalar (r.h.s.) form factor. 

multi-parameter fit and for the presentation here we discuss only the results of a linear fit. 

The derivative of the form factor at Q 2 = gives the charge radius, as shown in Fig. ||a. 
Whereas the large mass results are compatible with a constant, the lowest mass is smaller but with 
a very large statistical error. The values are quite compatible with values from at comparable 
quark mass values. A fit to a constant gives (r 2 ) v = 0.279(16) fm 2 . The current PDG average is 
0.45(1) g. 

The value of (r 2 ) v is in the VMD model inversely proportional to the mass squared of the 
p-meson. We obtain m 2 JT)M = 0.83(5) GeV 2 . This is close to the range of values obtained for m 2 
in the BGR analysis [§] in the direct p -channel (e.g., 0.80 GeV 2 at am = 0.04). These large values 
for the m p mass might at least in part explain (via the VMD model) our low result for (r 2 ) v . 

The scalar form factor also has disconnected contributions which we disregard here (like 
it is done usually due to the related technical complications). We determine the scalar radius 
squared from a linear fit to the <2 2 -dependence of the scalar form factor (now explicitely normal- 
ized according to Eq. 1.4). The resulting values, shown in Fig. |2|b, are almost an order of mag- 
nitude smaller than for the electromagnetic case. We obtain an average for the scalar radius of 
(r 2 ) s = 0.034(6) fm 2 . It is also smaller than the values expected for full QCD. 

In chiral perturbation theory the scalar radius is related to the quark mass dependence of the 
pion decay constant 

1 1 3 M 2 

hlf = 1 + ^ (r 2 )M + + 0{M A n ) , (4. 1) 

and in Ref. [||] values around (r 2 ) s ~ 0.55 — 0.75 fm were expected. Results for the pion decay 
constant in a recent full QCD lattice calculation []8|] lead to the value (r 2 ) s = 0.5 ±0.1 fm 2 . A 
corresponding analysis of the quenched BGR data [Q] gives again a small value 0.08 — 0.13 fm 2 . 
This seems to imply that (r 2 )^ is very sensitive to quenching. 



5. Conclusions 

Chirally improved fermions provide a valid framework for a first principles lattice study of the 
hadron structure. They are Ginsparg-Wilson-like fermions and less expensive than other types of 
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chiral fermions. We have begun the study of hadron structure by investigating the form factors of 
the pion. The vector form factor is consistent with expectations in its general form, although the 
charge radius is 40% smaller than in experiment. For the scalar form factor the scalar radius is 
much lower than one would expect from unquenched calculations. Both effects may be signals of 
quenching. 

Further progress could be achieved - as usual - by using better statistics, more momenta, and 
larger lattices, both in lattice units (to study finite volume effects) and in physical units (to access 
lower transferred momenta). 

Most computations were done on the Hitachi SR8000-F1 at the Leibniz Rechenzentrum in 
Munich, and we thank the staff for support. 
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